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ABSTRACT 


In  this  paper  a hybrid  Remes-differential  correction  algorithm  for  computing 
best  uniform  rational  approximants  on  a compact  subset  of  the  real  line  is 
developed.  This  algorithm  differs  from  the  classical  multiple  exchange  Remes 
algorithm  in  two  crucial  aspects.  First  of  all,  the  solving  of  a nonlinear 
■system  to  find  a best  approximation  on  a given  reference  set  in  each  iteration 
of  the  Remes  algorithm  is  replaced  with  the  differential  correction  algorithm 
to  compute  the  desired  best  approximation  on  the  reference  set.  Secondly, 
the  exchange  procedure  itself  has  been  modified  to  eliminate  the  possibility 
of  cycling  that  can  occur  in  the  usual  exchange  procedure.  This  second 
modification  is  necessary  to  guarantee  the  convergence  of  this  algorithm  on 
a finite  set  without  the  usual  normal  and  sufficiently  dense  assumptions 


I 


1.  Introduction 

This  paper  is  divided  into  two  parts.  In  the  first  part  we  consider  X a 
compact  subset  of  the  real  line  with  card(X)  >_  n + 2.  Let  C(X)  denote  the 
class  of  all  continuous  real  valued  functions  defined  on  X,  normed  with  the 
uniform  norm,  i.e.,  ||f||  = max{|f(x)|  : x £.  X}.  Let  n be  a positive  integer 

and  set 

R^(X)  = {r  = 1/p  : p £ Iln,  p(x)  > 0 for  all  x e X} 

where  11^  denotes  the  set  of  all  algebraic  polynomials  of  degree  ^_n.  Note 
that  R^(X)  consists  of  only  the  positive  elements  of  the  set  usually  denoted 
by  R^(X).  In  this  setting  we  will  give  an  algorithm  for  computing  the  best 
approximation  fcr  positive  f £ C(X)  from  R^(X),  We  believe  this  algorithm  is 
the  correct  analog,  for  this  setting,  of  the  standard  multiple  exchange  Remes 
algorithm  for  polynomials.  We  observe  here  that  if  YC  X,  Y is  compact  and 
card(Y)  n + 2,  then  existence  of  a best  approximant  to  positive  f €.  C(Y) 
from  R®(Y)  is  guaranteed  by  [5]. 

This  algorithm  contains  some  unique  features  including  the  incorporation 
of  the  differential  correction  algorithm  [1],  [4]  to  obtain  a best  approximation 
at  each  stage.  This  insures  that  the  denominator  of  the  best  approximation,  p , 

K 

on  the  ktn  reference  set,  X , will  be  positive  on  X,  . If,  however,  p,  (x)  < 0 

k k k — 

for  some  x £ X * X^,  we  indicate  two  exchange  procedures  for  selecting  the  next 

reference  set.  Note  that  in  most  studies  this  possibility  is  ignored  by 

assuming  (1)  f is  normal  on  some  interval  [a,  b]  containing  X;  (2)  X is  sufficient 1 

dense  in  [a,  b];  and,  (3)  X,  is  sufficiently  close  to  an  alternating  set  of  the 

best  approximation  to  f on  X.  We  shall  also  show  that  using  our  exchange  procedure 

there  exists  a k.  > 0 such  that  for  k > k , p,  must  be  positive  on  X.  From  this 
U 0 ^ 

point  on,  our  exchange  procedure  will  coincide  with  the  standard  multiple  exchange 
procedure  and  wc  can  therefore  guarantee  convergence  without  the  above  assumptions. 
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Our  procedure  could  also  be  used  to  overcome  the  difficulty  which  Dunham  [3] 
has  pointed  out  in  William's  paper  on  interpolating  rationals  [7]. 

It  should  be  further  emphasized  that  a modified  exchange  procedure  is 
actually  necessary  to  guarantee  the  convergence  of  this  algoi’ithm  without 
the  assumptions  ( 1 ) — ( 3 ) of  above.  Indeed,  if  one  attempts  to  use  the  standard 
exchange  procedure  without  regard  to  the  possibility  that  p.  < 0 on  X • X may 

K K 

occur  (and  hoping  that  p (x)  = 0 for  x € X ~ X does  not  occur  to  give  a 

K K 

divide  fault)  the  usual  proof  that  the  error  of  approximation  on  the  successive 
reference  sets  is  strictly  increasing  is  false.  In  fact,  examples  exist  for  which 
the  error  does  not  increase  strictly  and  for  which  the  algorithm  actually 
cycles  (i.e.,  ; pkt5  = = . . . ; Pkrl  = Pkt3  = Pk.5  = = \,2  = \t4 

= ...  ; = = = s'tartin2  a't  some  k).  Using  either  of  the 

exchange  procedures  that  we  give,  we  are  able  to  prove  that  the  error  of 
approximation  or.  successive  reference  sets  is  strictly  increasing. 

The  second  part  of  the  paper  is  devoted  to  the  description  of  the  Remes- 
Difcor  algorithm,  (the  name  of  our  algorithm)  for  obtaining  the  best  approximation 
to  f 6 C(X),  X a finite  subset  of  the  real  line  of  at  least  n + m + 2 points, 
by  elements  of  R™(X),  m _>  0,  n _>  0,  where 

R™(X)  = {r  = p/q  : p c nm,  q e TI^ , q(x)  _>  e for  all  x e X), 

and  the  e.  is  chosen  so  that  a best  approximation  from  R™(X)  will  also  be  a best 
approximation  from  the  larger  class  that  relaxes  this  requirement  to  q(x)  > 0 on 
X.  A proof  of  the  convergence  of  this  algorithm  is  given,  along  with  a flov; 
chart.  Finally,  a brief  discussion  of  some  numerical  results  will  be  given,  A 
complete  discussion  of  the  numerical  results  and  comparison  with  both  the  Remes 
algorithm  and  the  differential  correction  algorithm  is  planned  for  in  a separate 


paper. 
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2 . Approx. mat i r.g  t: 1 1 h R ( X ) . 

Let  f £ C(X)  ~ R°(X),  with  f > 0 on  X.  We  first  consider  the  case  where 
n 

X is  a finite  subset  of  the  real  line,  with  card  (X)  ^ n + 2.  For  each  k, 

k = 1,  2,  ...»  X,  c X shall  denote  a reference  set  of  n + 2 or  n t 3 points  and 

0 0 
r.  = 1/p  £ R (X  ) will  denote  the  best  approximation  to  f on  X from  R (X,  ). 

k k n k k n k 

This  best  approximation,  r,  , is  obtained  by  using  the  differential  correction 
algorithm  applied  to  the  point  set  . There  are  three  advantages  to  finding 

r via  the  differential  correction  algorithm  rather  than  via  solving  a nonlinear 

k 

system  of  equations:  a solution  is  guaranteed,  we  are  assured  that  > 0 on 

, and  no  extra  complications  will  arise  if  has  n + 3 points.  After  computin 
r , if  has  n + 3 points  we  delete  one  point  of  X^  to  get  a new  set  of  n + 2 

points,  taking  care  that  f - r,  alternates  on  Y . If  X,  consists  of  n + 2 

points,  then  we  set  = X^ . Set  e^  = max{|f(x)  - r^(x)|  : x£  Y.^},  Z^ 

= (x£X:  p,^(x)  > and  consider  the  following  two  exchange  procedures 


i . 

I 


for  constructing  the  next  reference  set  X. 


k+1  * 


k! 


t 


ip 

k 


Exchange  I : (The  positive  exchange)  If  r^  is  not  the  best  approximation  to  f 

on  Z,  from  R°(Z,  ),  X,  is  constructed  from  Y by  doing  an  ordinary  Remes 
K n k K + X k- 

multiple  exchange  on  the  points  of  Z^ . If  r^  is  the  best  approximation  to  f on 
Z^  from  R^CZj^.)  then  the  algorithm  terminates  if  Z^  = X.  If  Z^  i X then  y £ X 
satisfying  p^.(y)  = minfp^x)  : x £ X)  is  found  and  X^+-L  is  defined  to  be 
5^,+l  = (y}U  { Y^.}  . Note  that  in  this  case  we  have  the  p^ty)  0 and  X^+1 
consists  of  n + 3 distinct  points  of  X. 

Exchange  II.  (The  negative  exchange)  In  this  exchange  procedure,  the  algorithm 

first  does  a standard  Remes  multiple  exchange  on  the  point  set  Z^  with  respect 

to  f - r,  and  Y,  getting  W C Z.  , where  V.’,  consists  of  n + 2 points  on  which 

f - r^  alternates  in  sign,  |f(w)  - r^(w)[  >_  e^  for  all  w£  W^  and 

max{lf(w)  - r,  (w)|  : w £ W,  ) = naxf|f(x)  - r,  (x)|  : x £ Z.  ) . If  K = Y.  and 
1 f 1 k 1 k 1 x k k 


Z,  = Y and  Z,  = X then  the  algorithm  terminates  as  r is  the  desired  best 
k k k K 


approximation  to  f on  X.  If  this  does  not  happen  then  X is  defined  to  be 

K * X 

V<'k  U (y ) if  i X where  y satisfies  p^(y)  = min{pj<(x)  : x £ X}  0 and  W,  if 


Note  that  this  exchange  procedure  differs  from  the  first  one  in  that 

whenever  Z i X an  additional  point  were  p,  takes  on  its  minimum  is  added  to 
k k 

the  reference  set.  In  the  first  exchange  procedure  this  additional  point  is 
added  only  when  r^  is  the  best  approximation  to  f on  Z^  from  R^(Z^).  Also, 
note  that  whenever  Z^  = X both  of  these  procedures  coincide  with  the  standard 
Remes  multiple  exchange  procedure.  For  both  of  these  exchange  procedures  the 

following  theorem  holds.  (The  set  X^C  X is  chosen  so  that  it  has  n + 2 points 
and  e^  > 0 . ) 

THEOREM  1.  If  X is  finite  and  the  algorithm  described  above  using  either  of 
the  two  exchange  procedures  is  applied,  then  {e^}  is  strictly  increasing. 
Furthermore,  the  algorithm  eventually  terminates  at  a best  approximation  to  f 


on  X from  RU(X). 

n 

Proof:  To  show  that  e,^  < e^+^  for  all  k one  must  consider  two  cases.  The 

first  is  when  X^+^  is  constructed  only  from  points  of  Z^.  In  this  case  p^  and 

Pk+1  are  both  positive  on  Xj^  and  a standard  de  La  Vallee  Poussin  type  of 

argument  (zero  counting)  shows  that  e,  < e,  , since  p.  , is  best  on  X,  , and 
b h k k+1  k+1  M-l 

^ ^k+1’  ca3°  ^+1  = ^ {y}  where  = Y^  or  is  the  result  .of 

a standard  Remes  multiple  exchange  on  the  points  Z with  respect  to  f - r and 

k K 

Y , and  y £ X satisfies  p (y)  = minfp,  (x)  : x £ X)  <_  0,  there  are  two  subcases  to 
k k k 

be  considered.  The  first  is  when  !f(y)  - r.  ,(y)  < e,  , so  that  f - r,  , 

1 J k+1  k+1  k+1 

alternates  on  W,  with  error  e,  , . Since  p,  is  also  positive  on  W,  and 
k k+1  k k 

|f  - p | > e on  W we  must  have  that  e . ">  min(|f(z)  - r (z)|  : z € W ) e 
kkk  k'X  k k K 


by  the  same  de  La  Vallee  Poussin  type  of  argument.  Finally,  if  |f(y)  - r,  (y)| 

= and  f - r^+^  alternates  on  + 1 C.  where  Y^+^  t W^,  then  we  must 

have  that  max{]f(z)  - r^+  (z)|  : z £ } = e^+, . Also,  f - r^  alternates  in 

sign  on  W with  If  - r I > e,  on  W . Thus,  by  zero  counting  we  must  once  again 

k k k K 

have  that  max{ |f(z)  -r,  (z)  | : z 6.  W } >min{|f(z)  - r (z) | : z G W } since 

K-t-  J.  k K k 

p £ p implying  that  e,-  < e . (For  a more  careful  treatment  of  the  de  La 
k K + l K K r .L 

Vallee  Poussin  type  of  argument  see  the  proof  of  Lemma  2 later  in  the  paper.) 

The  rest  of  the  theorem  now  follows  since  X is  finite,  and  no  reference  set 
can  occur  more  than  once. 

Although,  in  actual  computation  one  only  encounters  finite  sets,  it  is  of 
interest  to  consider  the  behavior  of  this  algorithm  if  X is  only  required 
to  be  compact.  In  the  remainder  of  this  section  we  shall  only  consider  Exchange 
I (the  positive  exchange).  It  can  be  shown  that  similar  results  are  true  for 
Exchange  II.  We  first  note  that  in  this  case  the  set  = {x  <E  X : > O 

may'  fail  to  be  compact.  If  this  happens  then  it  may  not  be  possible  to  carry 
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no  g-poles  on  X for  k >_  some  ky  and  converge  uniformly  to  the  best  approximation 
r*  to  f on  X according  to  an  inequality  of  the  form  j|r^  - r:':||  ^ <_  A0^,  0 < 0 < 1, 
for  k ky 

Proof : Since  the  conclusion  follows  trivially  if  the  algorithm  terminates,  we 

oo 

assume  that  this  is  not  the  case.  The  method  of  proof  is  to  show  that  {e^},  , 

X K = 1 

is  increasing  and  to  actually  estimate  this  rate  of  increase.  To  prove  that 
e^  < ej<fq  holds  for  all  k,  one  simply  uses  the  arguments  of  Theorem  1,  Also, 

OO 

note  that  {e,  },  is  bounded  (otherwise  r = 1 would  be  a better  approximation 
than  rj,  on  for  some  k).  Hence,  there  exists  e*  such  that  ej,.  t e*.  The 
remainder  of  this  proof  is  broken  into  seven  lemmas;  the  first  of  these,  which 
proves  that  the  points  in  cannot  cluster  is  proved  by  arguments  similar  to 
Wendroff  [6,  p.  65]. 


k k 

LEMMA  1 . There  exists  6 > 0 such  that  for  every  k,  if  Y,^  = (Xy,  ...»  x +.,}, 


then  x*  <_  x^+^  - 6 for  i = 0,  1,  . ..,  n - 1. 


V 

l 


<r* 

K 


Proof : Suppose  not,  then  there  exist,  for  some  fixed  i,  subsequences  (relabelled 

\ k k k . k 

as)  {x.}  and  fx'.  ,}  such  that  x.  -+  xv  and  x.  , -+■  xv.  By  passing  to  further 
i l+l  i i l+l  i J v b 

V 

subsequences  (relabelling  if  necessary)  we  have  that  x'^  ->  xv  for  j=0,...,n+L 

as  k + co  where  xv  = x*  , and  x v < xv  , , i = 0 , 1 , ....  n - 1 . Thus , on  the 

l l+l  3 — 3+1  J 5 

set  X*  = ( x;'f , ...,  X"  , ) we  can  find  p*  £ n such  that  f(xv)  - — 7rr~^r  = 0, 

0 n+1  n 3 p-(xv) 

3=0,  1,  ...,n+l.  By  continuity,  there  exists  a 6 > 0 such  that 
e2  n+1 

|f(x)  ” ~yr( — t|  < n—  for  x £ U (xv  - 6,  x*  + 5)  P,  X where  > 0,  is 

* j=o  3 3 

the  error  of  the  second  cycle.  Hence  for  sufficiently  large  k,  we  have  that 
1 e ^ e,; 

If(x^)  - | j = 0,  1,  . ..,  n + 1.  But  this  implies  that  e.  < -37-  since 

1 1 p-(x!> ) 1 2 J k — 2 

1 . 3 j,-  k 

— is  best  on  Y^  = (xy,  ...,  x^-P  which  contradicts  the  fact  that  e^t.  l 
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LEMMA  2.  Let  X be  a compact  set  of  real  numbers  containing  at  least  m + n + 2 

points,  and  let  f £ C(X).  Suppose  r*  = £R™(x)  has  defect  d = nin(m-3p;'!,  n-8q:':) 

N — 

and  let  N=m+n+2-d.  Suppose  that  f - r*  alternates  in  sign  on  C X 

where  x^  < x^  < ...  < x^,  and  that  f(x^)  - r:':(x^)  i 0,  for  i = 1,  . ..,N.  Then 

if  r = — £ Rm(X),  r f r:’:  on  X we  have 
q n 

max  |f(x.)  - r(x.)|  > min  |f(x.)  - r:’:(x,  ) | . 
l<i<N  1 1 l<i<N  1 1 

Proof:  Suppose  max  lf(x.)  - r(x.)|  < min  |f(x.)  - r*  (x.)  . Let 
l<i<N'  1 1 ~ l<i<N  1 1 

A(x)  = r(x)  - r*(x)  = (f(x)  - r'-(x))  - (f(x)  - r(x)),  for  all  x £ X.  Assume 
(without  loss  of  generality)  that  f(x^)  - r*(x^)  > 0;  then  we  have  (-l)'A(x. ) <^C, 
i = 1,  . . . , N.  Mow  for  all  x £'  Y, 


""  J 'J  1 niin'W 
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6 

with  0 < 3^  < Let  J = I - (J  1^  and  let  6^  = m^n(  jq(x) | : x £ J}  >0, 


i=l 


62  ^ 

Choose  k so  large  that  |q  (x)|  > — for  all  x £ J and  c,  < 7-6  m where 

k — 1 — z z 

m = min{f(x):  x £ X}.  By  Lemma  1,  no  two  points  in  Y,  lie  in  the  same  I.; 


furthermore,  for  all  x e J,  we  have 

I — 7 — r|  = | — y-yj  < ^ m-^-  = m. 
pk(x)  - 2 2 o2 


, ! 


k k k 

Now  let  x.  , , x.  and  x.,,  be  consecutive  points  of  Y,  and  suppose  that 
1-1  1 1+1  k ™ 

k 1 k lk 

f (x • ,)  - ; = e,  . Then  x^  must  lie  in  some  I.  (since  ; — = f (x.)  + e,  >n) 

1-1  / k s k 1 i / 1 k 

P Xi-]  k k pk  Xi 

which  is  separated  from  both  x-  , and  x.  . by  ooints  of  J. 

1-1  1+1  J ■ 

Since  | — y- y|  < — at  such  separation  points,  if  p >0  throughout 

PkU;  Pk(Xj)  k 

b V _ m t 1C 

[x.  , , x.  ],  then  p must  have  a relative  minimum  somewhere  in  (x.  , , x.  ). 

• 1-1  1+1  k ]-l  i+l 

k k 

If,  on  the  other  hand,  p (x)  < 0 for  some  x £ [x.  , x.  ] then,  since 

k 1 - 1 1 + 1 

k k • 

p (x.  ,)  > 0 and  p, (x.  ,)  > 0,  it  again  follows  that  p,  has  a relative  minimum 
k 1-1  k l+l  k 

k k 

somewhere  in  (x.  , ,x.  ,)• 

1-1  1+1 

k k k . . 1 

Next,  assume  x.  , , x.  and  x.,,  are  consecutive  noints  of  Y,  with  

V1’  3 VI  - k <xk  ) 

k J-. 


f (xk  ) = e, ; then 
1-1  k’ 


p (x.  ) 

Fk  l-l 


= f(xk  , ) + e,  > m implying  that  xk  £ I.  for 

k •,  1-1  k 1-1  i 


y 


• • • k v k 

some  i,  and  similarly  for  xt  , . By  Lemma  1,  x.  , and  x.  , are  in  distinct  I.’s. 

1+1  1-1  1+1  l 

k k k k 

If  p^_  ^ 0 throughout  either  [x_.  ^ , x_.  ] or  [x_. , x^+^],  then  for  some  point  x in 


one  of  these  intervals  we  have  0 < -—7—*-  < — , and  0 < — t — r < 

P,  (x)  /,  k v p,  U)  / k 

k p k X j - 1 k Pk(xj  + 1} 

• . k k v k 

since  there  are  points  of  J between  x.  , and  x.  and  also  between  x.‘  and  >: . , . 

1-1  1 1 1+1 


inrr  ir>  r 
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k k 

Therefore,  p,  must  have  a relative  maximum  somewhere  in  (x.  ,,  x If,  on 
rk  3-1  3+1 

k k k k 

the  other  hand,  p.  < 0 somewhere  in  both  [x.  , x . ] and  [x.,  x.  ],  then  sir.ee 

K j-l  J 1 1+1 

k 

p (x.)  > 0 it  again  follows  that  p,  has  a relative  maximum  somewhere  in 
k 3 k 

, k k . 

(xj-l’  Xj+1}- 

We  have  now  shown  that  n,  has  a relative  minimum  between  every  pair  of 

- k 

"lower  extrema"  of  f (on  Y ) and  a relative  maximum  between  every  pair 

pk  k 

of  "upper  extrema".  Thus,  pk  has  at  least  n relative  extrema.  But  p^  is  a 
non-trivial  polynomial  of  degree  <_  n.  This  contradiction  completes  the  proof 
of  the  lemma . 


Corollary.  There  exists  a constant  c:':  > 0 such  that  |pk(x)  [ <_  c*  for 
k = 1,  2,  ...  and  all  x 6.  X. 

Before  proceeding  to  Lemma  4,  we  introduce  some  new  notation  and  make  a 
few  remarks.  We  shall  call  the  exchange  from  to  X^  an  augmented  exchange 
if  X,  = Y,  U (y)  (recall  that  Y,„  C X,  is  a set  of  n + 2 points  on  which  f - r, 

KM  K K K K 

alternates  with  error  e,  ).  Also,  in  this  case  the  point  y £ X satisfies 

k k 

Pk(y)  = min{pk(x):  x £ X}  < e.  Writing  Xk+1  = {yQ , ...,  Yn+2 K we  have  that 

k 

X,  , contains  exactly  one  g-pole  of  r.  . Call  this  point  y . As  stated  earlier, 
k+1  K a 

we  let  r,  , denote  the  best  approximation  to  f from  R^(X,  , ) on  X,  , (found  via 
k+1  n k+1  k+1 

the  differential  correction  algorithm)  and  we  define  Yk  to  be  that  subset  of 

X,  , on  v;hich  f - r,  , alternates  in  sign  with  modulus  e.  Note  that  since 

k+1  k+1  k+l 

we  are  assuming  that  wc  are  using  Exchange  I,  Y is  uniquely  determined  by 

k 1 1 

the  fact  that  yQ  c Y ^ and  at  precisely  one  point  of  Yk,  say  t,  |f(t)  - rk+^(t)l 

must  hold.  This  follows  from  Lemma  2.  For  e.  < A < e,  , construct  r,  = 

k — — k+1  A 

by  requiring  that 

f(y^)  - rx(yk)  = n^X,  i = 0,  1,  ....  n 


&& £***&&£ 
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where  Xk+1  = - (y^>  = <V  } 3nd  ni  = 

~k  -k 

Observe  that  for  = -1,  f(y^)  - ruA  = f(y^)  + A > 0 and  for  = +1, 

f(y^)  - = f(y^)  - x 1 f(y^)  - ek+1  = rk+1(y^)>  \ > 0.  Thus,  pA  is  wei: 

defined  by  these  equations  and  for  all  x G.  X,  p^Cx)  is  a continuous  function  of 
Finally,  let  A = inf{  ||f  - r||  : r £ K^(X)}.  Note  that  A > 0 since  f C C(X) 

- R°(X).  Then, 


LEMMA  4 . If  at  the  k-th  exchange  an  augmented  exchange  occurs  and 

sgn[f(yk)  - rk(yh]  = sgn[f(y^)  - rk+1(y^>  i = 0,  1,  ...,  n,  then 

e,  - e,  _>  ft(  ||  f ||  - A),  where  ft  is  a constant  independent  of  k. 

Proof:  First  observe  that  if  X = e,  , then  r = rv  as  these  two  functions 


take  on  the  same  values  on  X,  (n  + 1 points)  and  likewise  r = r.  . 

KtI  6.  . KtJ. 

- k k 

Thus,  p (y  ) < e since  y is  a g-nole  of  r, . 

a J a & - k 

k , i ■, 


We  now  claim  th 


. c -k  k -k 
lat  if  y.  < v < y. 


y.  < v < y.  then  we  must  have  that 
i ' o i+I 


sgn[f (y ^ ) - rktl(yh]  = 3gn[f(yJ+1)  - ^(y^)]  = 1 


since  otherwise  (i.e.,  = -1)  a zero  counting  argument  implies  that  r,  = r.  . 

K KtJ. 

k -k  -k  ■“  k ^k  ~k 

Similarly,  if  y'Q  < yQ  or  y^  < y , we  must  have  sgn[f(yQ)  - rk+1(yQ)]  = 1 or 

'k  ~k 

sgn[f(y“)  - rka(yn)]  = 1,  respectively.  This  follows  by  counting  zeros  of 

Pk+j  - pk  in  the  contrary  case  and  using  the  fact  that  ek  < ek  . Indeed, 

-k  v 

suppose  y < y"  and  sgn[f(y ')  - r,  , (y  )]  = -1.  Now  vre  have  that  r,  , alternates 
n Jo  b Jn  k+l  Jn  k+1 

~k  ~*k  k k ^ k 

on  y',  ...,  y ',  y„  with  error  e,  , and  r,  alternates  on  y'\  ...,  y with  error 
0 n 0 k+l  k 0 n 

e,  and  the  same  sign  as  that  of  f - r,  . Thus,  d,  , - p,  has  n zeros  in 
k k+l  - k+l  rk 

^ k ~ k ~*k  ~V 

[y«,  y‘].  Also,  we  must  have  r,  (y;  ) > r (y  ) since  we  are  assuming  that 
v n k i"x  n k n 

~k  ^k  -k  -k 

sgn[f(yi)  - rk+1(yi)]  = -1  and  ek  < ek+1-  Thus,  Pkfl(yn'>  < pk(yn>.  Also,  since 

J, 

rk+1  is  the  bust  approximation  to  f on  Xk+1  we  must  have  Pk(.}(yJ)  > e • Now, 


irar 
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k k 

by  construction,  y was  chosen  so  that  o,  (y  ) < e.  Thus,  p.  - v has 

o - k a — rK  + l k 

~k  k 

(at  least)  one  additional  uero  in  (y  , y ) implying  p,  , = p,  which  is  a 

J n ’ a 1 J u ■ k+i  k 

contradiction.  A similar  argument  will  treat  the  other  cases.  Thus,  it 

k k 

follows  from  the  alternation  of  f - r,  , , that  sgn[f(y  ) - r,  ..(y';]  = -1. 

k + 1 o k tl  o 

k - k k 

Consider  p,(y  ) for  e,  < X < e,  , , . We  know  that  p (y  ) = p. (y  ) < e. 
*K  'a  k — — k+1  ' e,  Ja  k a 

Set  m = min{r(x):  x X}  and  ” = max{f(x):  x 6.  X}.  Since  |Jf|j  > 2A  > A 


we  have  (as  e < 


rimi 


1 ) 


■k+1' 


k — 

(y  ) = p / k.  k,  .-l 

0 ek+l(ya}  = (f(y«}  + ) 


k+1' 


> (f(y^)  + ek+1  + ( |jf||  - A))  1 


” (f(ya}  + ek  + l + ( HfH  " ek+l))-1 


i -nmi  i = ^ V#- 


Therefore,  p (y^><(f(yg)  + e + ( |jf]|  - A)  )-X  < p,  (v£ ) . Since  p is  a 


-1 


k c 


K + 


k+1  - o 


continuous  function  of  X,  e,  < X < e,  , , there  exists  w such  that  e.  < a,  < «>.  . 

k — — k+1’  k k + j. 

and 


i,o,(ya)  = (f(yo)  + \+i  + ( HfH  " &)) 


-1 


Wc  define  coefficients  c.  for  j = 0,  1,  . . . , n by  setting 

1 5 A 

PA(yi)  = l CA  X(yi)] 

A 1 j = 0 J> 

and  let  M*  = max{l,  |x|,  ...,  |x|n).  Then, 
xeX 


|jf||  - ^ = (ek^  + f(y*)  + ( ||f I!  - A))  - (ek+1  + f(y*)) 

k , 


PekH{yo)  - Pw(ya) 


k 


pJyV  p„  p Jy''P„ 


u ' o 


ck1i  a 


u o e.  ■ o 

kvd 


i-Up  (yak)  - p(i  (yk )] 

e2  ek+l  a 01  0 
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' 


- =i  >0ktb 

e ^ j=0  J,ek+1  J,u  0 


„ M*(n  + 1) 

< r max  c . 


- c 


, J 


e2  °£3in  ]®k+l  :,w 

M'-(ri  + 1)  i I 

x c . - c . 

e2  ?’ektl 


We  complete  the  proof  of  Lemma  4 by  showing  that  jjf||  - A <_  ft'(eR+^  - e^). 


•where  is  a 

constant  independent 

of  k . Now 

let 

k 1 1 9 ° 

, ~k,2 

(y0)  •• 

,~k.n 

(y0) 

D(y£)  = det  : • 

• 

• 

; 

\ : ~k 

V 1 yn 

,ik,2 
(yn;  •* 

,ik,n 

yn 

and  let  D(y\ 

k)  be  D(y^)  with  the 

j-th  column 

replaced 

((f ($q>  - V)_1’  (f (^n)  - nnX)’1)T 


Also,  let  K.  be  the  cofactor  of  D(y*,  M relative  to  the  (i,  j)  element, 
using  Cramer's  rule 


I hen 


«■> 

N 


Ilf  II  - 4 < : 


|D(y„k)l 


M* ( n t 1 ) 
k. 


I w.  C- 


iCD(ya’  ek+l)  ‘ D(V  u,)]! 
1 1 


r ~k, 


e2|b(yX)|  i=0  * f{yp  - nieR+,  f(y^)  - ru. 


1 


— 3 

U) 


< HMnJh  12  l |W.|[m,x r 

e2|D(yX)|  i=0  1 °-L-n  If(y^)  - n.ek+1l  |f(yj)  - I 


](fck+1"“) 


''k 

Now,  since  the  points  y^  are  separated  (Lemma  1),  there  exists  a constant  f,  * 0 

(independent  of  k)  such  that  the  Vandermonde  determinant  ! D ( y^ ) I f > 0. 

n 

Furthermore,  £ |w.|  K for  some  positive  constant  K since  all  cofactors  cf 
1 i=0  1 

D(y0  ) are  bounds]  (independent  of  k).  Finally, 


* 

* 


I 


!f<n>  - vkJ  ■ - hr,  - —7^7  i & 

V.  ,‘yi'  Pktl^i’ 

k+1 

by  the  corollary  following  Lemma  3.  Furthermore,  |f(y^)  - n-w|  = — — ; — , 

~k  1 p (y.) 

and  by  construction  this  lies  between  f(y • ) and  , so  that  u 1 

Pe  (yj) 

|f(y^)  - n • W ! > min(m,  -4r) • k+1 

1 j.  1 — e#* 

Thus, 


|f(yj)  - Vk+1|  |f(y*)  - v! 


< o*max(e'“,  — ) . 
— m 


Therefore , 


- A < ft*(e,  - u).  Taking  f2  = 777-  yields 


n(  llf|1  - A)  1 ek+i  " “ 1 ek+i  - V 


Now  let  us  assume  there  exists  a subsequence  of  positive  integers 
00 

{km},n_1  satisfying  the  following: 

1.  An  augmented  exchange  occurs  between  Xy  and 

"m  m+1 

2.  e^  +1  - < fj(  j|f  ||  - JJ)  (since  e^  + e*  _<  A,  where  A is  the 

m " m tn 

error  of  best  approximation  to  f from  ( X ) ) . 


By  our  assumption  2,  we  see  that  the  sign  condition  of  Lemma  4 cannot  hold, 

hence  we  have  for  each  k , the  additional  condition: 

m 


3.  sgn[f(y.m)  - r,  (y.m)]  i sgn[f(y.m)  - r (y.m)]  for  some  i,  0 <_  i < n. 

2.  K t J.  1 1 K 1 


Recall  that  for  an  augmented  exchange  between  X,  and  , that  denotes 

mm  'm 

the  subset  of  X,  consisting  of  n + 2 points  on  which  f - r.  alternates 
K t 1 K < J. 

m m 

k k k 

with  error  +1 . Define  yg  ‘ by  (y™)  = Yk  ~ Yy  +1-  That  is,  y™  is  the 
m m 'm 

point  of  a,  which  is  deleted  in  forming  Y,  , . Since  we  are  considering 

k +1  k+1 
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Exchange  procedure  I,  we  have  that  W may  be  taken  to  be  Y whenever  an 

k k 

m m 

augmented  exchange  is  performed.  Under  these  assumptions  we  prove  the 
following  two  lemmas . 


— k g 

LEMMA  5.  If  p (y  ) < 7- r-  , then  e,  , - e,  > S2"  where  ft"  is  a constant 

ek  ^ (1  + k+1  k ~ 

independent  of  k. 

' Proof:  Once  again,  define  r^  = X-  , e^  £ A <_  ek+^  by 

px(y^)  = (f(y^)  - niek)"1,  0 < i < n 
-k  -k 

with  n-  = sgn[f(y.)  - r,  ,(y-)]»  0 < i < n.  Since  we  are  no  longer  assuming 

1 1 Kt X X 

^k  *»k  ~k  -k 

that  sgn[f (y . ) - r,  (y - ) D = sgn[f(y.)  - r (y.)]  for  all  i,  we  do  not  necessarily 
1 Kl  1 Kt  JL  1 

have  that  p = p,  . However,  we  still  have  that  p = p as  before  so  that 
. ek  k ek+l  k+1 


Pek+I(yke)  I G > ufry  holds-  since  Pek(ye)  " tttit 


by  hypothesis  we  have. 


by  the  Intermediate  Value  Theorem,  that  there  exists  m,  ek  < w < ek+^,  such  that 

pto(ye}  = ("~t  'el  • Then’  as  before 


! = XJL£  . I < 

e s — 


- , k.  - , k. 

p„<ye>  >Vl<V 


?ek,i(ye>  - 


t' 


< [ (c.  - c.  Ky^ 

~e  j=0  ]’ektl  6 


where  the  coefficients  c. 


and  c . are  as  defined  in  the  proof  of  Lemma  4 . 


l,ek+i  D,co 

Since  the  estimates  used  in  the  proof  of  Lemma  4 are  independent  of  the  point 
k k 

y^,  they  may  be  applied  here  for  y^,  and  we  have  1 <_  ftg(ek+^  *'  <*>)  1.  ^Q^ek  + 1 ' c'k 
implying  ek+1  - ek  > j-  = ft"  > 0. 


e,  ) 


*7* 

4 

[ * 

La 
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LEMMA  6.  An  augmented  exchange  can  occur  only  a finite  number  of  times. 

Proof : Assume  the  contrary.  Then,  by  extracting  subsequences  (as  often  as 

k ^ ^2, 

necessary)  we  obtain  a sequence  {y^  } of  extreme  points  (xq,  ...»  xn+^)  such 
that  the  exchange  from  Y,  to  X.  is  an  augmented  exchange,  so  that  X,  . 

k£  k£+1  V1 

, k»  kj, 

= Y,  (J  (y  ) where  y is  a g-pole  for  r,  and  is  selected  so  that 

k«,  0 ° k2, 

p (y  ; = min{p,  (y):  y G X}  < e;  r.  is  the  best  approximation  to  f on 

k 0 

Z,  = (x  G X:  p,  (x)  > e}  from  R°(Z,  ) so  that  no  multiple  exchange  is  applied 

k£  k£  n k£ 

k2.  k£ 

to  Y,  . Letting  y be  defined  as  before,  {y  }=  Y - Y , we  assume  that  the 

k£  0 2 k£  V1 

final  subsequence  (Y  } for  k = 1,  2,  ...  satisfies 

k£ 


(1>  \ (y6  > i (TT7T 

k£ 


( Lemma  5 ) . 


(2)  sgn[f  (y^)  - rk(yiV')]  i sgntfty^)  - r,+1(y.1)]  for  some  i,  0 £ i £ n, 

^ * kj>  ko  k«  . , 

independent  of  k , where  X^  +1  = Y^  +1  - {y^  } = {yQ  , . . . , yR  } (Lemma  4;. 

I 2,  0 

kp 

(3)  y •+  y’£  G X as  k -+■  oo 

6 & £ 


(4)  X X*  = {y*  ...,  y *}  c X (coordinatewise  convergence)  with 

k tx  u n 

% 


y>5  , - y v > 6 > 0 , 0 < j < n ( Lemma  1 ) . 
J3+l  J - 


(5)  p,  -►  p-  €.  n uniformly  on  X (Lemma  2). 
rk  n 

l 


**  X,  j • • 

(6)  sgn[f(y^  ) - r^  (y..  ')]  is  constant  for  fixed  j,  independent  of  kj 


As  noted  above,  all  of  these  conditions  can  be  met  by  passing  to  subsequences  of 
subsequences  sufficiently  often.  Now,  under  ’’hese  conditions,  we  claim  there 


m y>«n 


13 


exists  a p > 0,  independent  of  k^,  such  that 


. k p _ . 

lf(y3  > - re,  (yp  )l-ek/p 


(2.1) 


where  r,  is  defined  for  e,  < X < e,  as  before: 

x k*-  - V1 

„ k g ^ k C 

r\(y.  ')  = f(y.  ) - n-A.  3 < j < n 
A J J J 

_k  -k^ 

with  n-  r sgn[f(y.  ) - r,  (y.  )].  Indeed,  if  (2.1)  is  not  true,  then  there 
1 ‘I  Y1  3 

exists  a subsequence  (relabelled)  such  that 


Define  q €.  Hn  by 


(2.2) 


k5  -1  - _kn 

q(yv)  = lim  (f(y .l)  - q.e,  ) = lim  p (y.fc),  0 < j < n.  (2.3) 

3 • *r  3 3 \ kr  \ 3 

k 0 „k  n 

Note  that  if  n . ~ -1,  then  f(y.  ) - >_  f(y.  ) >_  min{f(x):  x c X}  = tn  > 0, 

3 3 3 l 3 

and  if  nj  = 1,  then  f(y^£)  - nj=k  > ~ \ +1  = vk  +i(yj£)  >■  ^ > °‘ 

£ & Z 

Thus,  the  above  limit  exists  for  each  j and  q(x)  > 0 on  X*.  Prom  this  it 

follows  that  po  converges  uniformly  to  q on  X.  Furthermore,  q(y*)  >_ yr-f — c 

by  (1)  and  by  (2.3)  [f(y::)  - —lev I = e*.  Thus,  for  k sufficiently  large, 

3 qky^w  Jo 

kj_ 

say  k.  > k,  so  that  q(y  ) > 0 we  have 

x.  p 

l«y|>-5<ir)-liU(y|)-f<y^)|tifCy>-?  (y^)| 

f ifev  (v>  - 7^r;i + - 5T7P1  * s* as  v 

q(ye  ) q(y6  ) 3 

,k^  kf 

Now,  since  p,  (y.  ) > e for  all  j and  k and  p,  (y  ) > e for  all  k , we  have 
k,  ] — x.  k.  a — z 


that  p*(y*)  > e for  all  j and  p*(y*)  > e.  Furthermore , f - — v-  alternates  on 
1 ~ 3 — p-‘ 


19 


1 " k 9 

X*  U {y*}  with  deviation  e*  since  f - alternates  on  X U {y  J with 

P ^ ^ 

1 l l 

deviation  e,  . Thus,  by  Lemma  2,  — = — . 3ut  this  is  impossible  since  for 

J q P" 

k^,  sufficiently  large  (i  the  index  of  (2)) 


sgn[f (y‘v)  - — qy-y]  = sgn[f(y.  ) - r^y..  )] 

~k  o .ki  ~k  t _ .kn 

= -sgn[f (y . ) - r (y.  )]  = -sgn[f(y.  ) - r (y.  )] 

1 kri  1 1 ekA+i  1 

= -sgn[f (y| ) - ^y]  i 0, 

which  is  our  desired  contradiction.  Thus,  (2.1)  holds.  Since  e > e,  + e::, 

V1  \ 

we  have  that  e - e,  -*■  0.  Choose  k so  that  k > k implies  e,  + £■  > e . 

V1  k£  £ ~ i V1 

Then,  we  have,  us ing  (2.1)  and  the  determinant  argument  from  the  proof  of  Lemma  4, 


O ! ko  _ kn  ko  _ kn 

2 < ek„  + P - ek„  + l - lf(ye  } - re,  (y6  }|  “ |f(ys  } ' re,  <yS  5 


kj  - ko 

pe  (y6  > pe  (y3  > 

kA  V1 


k o _ k o 

>e  (y3  } ' Pe  (yS  M 

V1  8 % 8 I 

\(yB  ’ • V6  ’ 


1 + e - / n / \ 

V 8 \ 6 1 1 V1'ek‘ 


Thus,  e - e >_  ft  , ft„  = • But  this  is  impossible,  so  we  have  that  an 

K +1  K Z Z.  Z^L  - 

i i 1 

augmented  exchange  can  occur  only  a finite  number  of  times.  I 

We  now  turn  our  attention  to  the  case  that  the  exchange  from  Y,  to  X,  , is 

k k+1 


not  an  augmented  exchange.  In  this  case,  r^  is  not  the  best  approximation  to 

f on  Z,  from  R^(Z.  ) and  X,  , = Y,  = {xk+\  ...»  x”^}  with  no  g-pole  of  r in 

k n k kil  Kt  1 0 n+l  k 
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v . Setting  y = min  |f(xj+1)  - r,(xk't'1)|  and  8 = max  | f (xk  "Vr,  (xr‘  + i . 

Kti  kti  0<i<n+l  1 ^1  k+1  0<i<'n*tl  1 .<i 


we  observe  that  yk+1  1 «k  and  3k+1  > ek+r 


LEMMA  7.  There  exists  a constant  S2  > 0 (independent  of  k)  such  that  if  X,  . is 

Kt  j. 

not  obtained  by  an  augmented  exchange  then  e,  , - y,  > J2(  8,  , - e,  ) . 

k+.l  'k+l  — k+1  k+1 

Proof:  Let  X be  a parameter  satisfying  Yk+1  £ X £ ek  . n 

= sgn[f(xk+1)  - 1'k(xQ+'*’)]  and  note  that  (-l)^n  = sgn[f(xk+")  - rk(xk+')]  for 

k+]  k+1 

l = 0,  1,  n + 1,  In  addition,  it  is  always  true  that  sgp[f(x.  ) - r,(x.  )] 

1 K 1 

= sgn[f (xk+k ) - r (xk+1)].  This  fact  follows  from  a zero  counting  argument  since 
l k+1  l 

both  r,  and  r,  , , are  positive  on  Y,  , , and  both  f - r,  and  f - r,  , alternate  in 
k k+1  k+1  k k+1 

l k+1  1 k+1 

sign  on  •k+-L*  If  (-1)  n = 1,  then  we  have  that  f(x?  ) - (-1)  nX  = f(x^  ) - X 

k+1  k+1  1 

>_  i(x^  ) - ek+1  = rk+^^xi  ) d.  C*  > 0.  Cn  the  other  hand,  if  (-1)  n = -1, 

we  have  that  f(xk+k)  - (-l)^nX  = f(xk+k)  + X >_  f(xk+k)  >_  m = mi.n{f(x):  x a 0 > 0. 

In  either  case  we  have  that 

k + 1 i 

f(x.  ) - (-1)  nX  > 0 for  0 < i < n + 1,  y,  , < X < e,  , . 

l - k+1  - - k+1 

Define  p,  t II  by  p.(xk+1)  = (f(x'.'+1)  - (-l)1nX)  1 for  i = 0,  1,  . n + 1, 
a n a i i 

i t q,  where  q is  the  smallest  subscript,  0 <_  q < n + 1,  for  which 

|f(xk+1)  - r.  (xk+1)|  = 3,  , • Next,  define  r,  = 1 and  note  that  r = r,  , 

' q k q k+1  X - ek+1  ktf 

* k+1 

since  these  two  functions  agree  at  n + 1 points  (i.e.,  . i = 0,  1,  . n+1,  i^q), 

- k+ 1 

Finally,  observe  that  p^(>:  ) is  a continuous  function  of  A for  y.  , . < A < e.  , . 

A q k+1  — — k+ J 

- k+1 

We  shall  prove  that  there  exists  an  to,  y,  , < w < e,  , such  that  p (x  “)  = 

1 'k+1  — k+]  1 to  q 

P (xk+1).  To  do  this  we  must  consider  two  cases: 


-i  i (^rrrm*rnmnr-~  * r.  *-r 


Case  JL:  (-l)°‘n  = 1,  i.e.,  sgn[f(x^  ) - rk+1(xq  )3  = !•  Here 


mktl)  - rj<(xq+1)  = Bk+1  “d  ffXq’1)  - WsT'’  ■ Vr  ror  1 * q-  ”e 


have  |f(xk+1)  - r (xk+1)|  = Yk+1  £ |f(x^+1)  - rk(xk+1)|,  so  that 


(-l)Xn[r,  (xk+1)  - p (Xk+1)]  < 0,  and  thus  (-1 )xn[p.  (x^*1 ) - p (x^+1)  > 0. 
k 1 Yk+i  1 k 1 Yk+1  1 


i t q.  Now,  if  (-1  )%[pk( x*  X)  - p (x  )>  0 holds,  then  by  counting  zeros 


Yk+1  q 


(including  multiplicities  of  up  to  order  2)  one  has  that  p,  = p , so  that 

* Yk+1 


k+1  - k+1 

pk(x^  ) = p^  (x^  ) and  one  sets  u)  = in  this  case.  If,  on  the  other 

k+ 1 • 


hand  (-l)qn[p,  ( xk+1)  - p (x*+1)]  < 0 holds,  then  p (x*+1)  > p,(x*+1). 


Yk+i  q 


Yk+i  q 


k+1,  . k+1.  . k+1.  , k+1.  . 

Since  f(xq  ) - rk(xq  ) = 6^  > eR+1  = f(xq  ) - rk+]Uq  ) we  also  ha.e 


k+1  k+1  — k+1 

that  p,  (x  ) > p,  ,(x  ) = P (x  ) so  that  by  the  Intermediate  Value 

^k  q * k+1  q rek+1  q 


- k+1  k+1 

Theorem  there  is  an  w,  yk  < w < ek+j  such  that  p^(xq  ) = Pk(xq  )■ 


Case  2:  (-l)n  = -1.  This  case  follows  with  essentially  the  same  argument 


and  we  shall  not  give  the  details 


_ i k ’ 1 

Thus,  there  exists  an  u,  Yk+q  <_  u < ek+^’  such  that  p^(xq  ) = Pk(xq  )■ 


Hence 


k+1.  - , k+1.  r /- / k+1.  / k+1,.,  1 1 

3,  ,-e.  =f(x  )-r  (x  )-[f(x  )-r  (x  )]  = — r— t — 

k+1  k+1  q ai  q q ek+1  q - (x*+l)  - (^  + 1) 

6k+l  q W q 


= k y [c.  - c ](>k+1)^. 

- , k+1,-  , k+1,  / 1 j,w  q 

Po,(xq  pe,  Xq  > :=°  k+1 

Kt  X 


- k+1  k+1 

Since  p (x  ) - p,  (x  ) > e,  and  referring  to  the  already  established 
oj  q k q — 


estimates  of  Lemma  4,  we  have 


ft  - e < -i— — J |W.  |[max(c*,  -)](e  - oj) 

k+1  ktl  - e2  |D(  k+l}  | /q1  i * k+1 


I 
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Now  we  may  apply  the  argument  given  in  the  continuous  case  ([2],  p.  99), 
noting  that  we  have  a Strong  Uniquenss  Theorem  ([5,  Theorem  3]),  to  show  that 
there  exists  0 6 (0,  1)  and  A > 0 such  that  if  k > k^  then 
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3.  Approximation  from  Rm(X) 

— n 

We  now  turn  to  the  second  objective  of  this  paper.  Here  our  approximat ng 
family  is  taken  to  be 


R™(X)  = {r  = p c IIn,  q g.  IT^,  q > 0 on  X}, 


and  require  card(X)  >_m  + n + 2 Cm  >_  0,  n >_  0).  g-poles  are  defined  as  before, 

i.e.,  x £ X is  said  to  be  a g-pole  of  r = ^ if  q(x)  < e where  e > 0. 

This  concept  is  useful  even  when  X is  finite, 

since  it  enables  us  to  avoid  division  by  very  small  positive  numbers.  We 

“18 

have  used  e = 10  on  a UN I VAC  1106,  which  has  roughly  18-digit  accuracy  in 


double  precision.  Unfortunately,  we  can  no  longer  be  sure  that  r^  will  be  g-pole 
free  or.  its  reference  set,  although  this  condition  can  be  enforced  by  inser' ing 
additional  constraints  into  the  linear  programming  part  of  the  differential 
correction  algorithm  (we  will  return  to  this  point  later).  The  algorithm  \u. 


r 


4 • Convergence  of  t h o Re me s - D : f c: o ; 1 alr-rithn 

In  this  section  we  prove  that  if  the  20-st  ,-p  stopping  criterion  is 

deleted  from.  the  Remes-Dif cor  flow  chart,  then  under  certain  existence 

assumptions  the  algorithm  will  terminate  at  a best  approximation  to  f from 

Rm(X). 

n 


THEOREM  ? . Let  X be  a finite  set  of  real  numbers  containing  at  least 
m + n + 3 points,  and  let  f m C(X).  Suppose  that  for  each  subset  Y C X 
containing  exactly  m+n+  2 or  m + n + 3 poincs,  a best  approximation 

EH!  fn 

£ R (Y)  exists  for  f from  R''(Y)  and,  in  addition,  that  o > e on  Y. 

q n n — 

Then  the  Remes-Difcor  algorithm  will  terminate  at  a best  approximation  r:':  to 

f on  X from  Rm  ( X ) . 

n 

Proof:  Let  XQ  be  the  initial  reference  set  and  let  X^  be  the  reference  set 

at  the  k-t’n  stage.  Let  r^  be  the  best  approximation  to  f on  X^  with 

ej,  = max(|f(x)  - r^(x)(:  x£  X^}.  If  the  algorithm  terminates  at  stage  k, 

then  there  are  no  g-poles  and  the  maximum  error  occurs  in  X^;  thus, 

e.  = Ilf  - r,  ||  and  r,  is  the  best  approximation  to  f on  X from  R^CX). 

Now  suppose  the  algorithm  does  not  terminate  at  the  k-th  stage  (k  >_  1). 

If  has  g-poles  in  X,  then  at  least  one  of  these  is  included  in  X^  by 

construction,  so  that  r,  t r,  Also,  if  r,  . has  no  g-poles  in  X,  then 

k k— 1 k-„ 

r,  tv.  , ; since  otherwise  the  maximum  error  for  r,  in  X would  occur  at  some 
K k-tl  K 

point,  as  the  maximum  error  for  r^  ^ , and  thus  would  be  included  in  Xj..  This 

would  contradict  the  fact  that  the  algorithm  does  not  terminate  at  the  k-th  stage. 

Now  f - r,  , must  alternate  on  some  set  {x, , x ...,  x . , } C X, 

k-1  1 2 - m+n+2-d,  , k-1 

k-1 

where  d,  , is  the  defect  of  r,  , . and  so  by  construction  f - r,  , must  alternate 
k-1  k-1'  k-1 

in  sign  on  some  set  {x, , x’ , . . . , x’  , } C X,  with  |f(x!)  - r,  , (x!)|  > 0,  . 

0 1 1 m+n+2-d,  , k 1 i k-1  1 1 — 

k-1 
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5.  Examples  and  C-p.cIusions 

In  order  to  get  a time  comparison  of  the  Remes-Difcor  algorithm  with  the 
ordinary  differential  correction  algorithm  alone,  we  ran  the  following  digital 
filter  design  problem: 

I 

1,  0 < x < 0.?n 


Let  X = [0,  0.2ir]  'J  [0.4m,  it],  f(x)  = 


0.0123, 


0 . 4ir  < x < it 


p ( X ) 

We  approximate  from  R„(X)  = = (a  + a cos  x + ...  + anccs  9x)/(t>  + b,cos  x 

2 q(x)  01  9 01 

+ b^cos  2x)  : q > 0 on  X).  We  also  want  q > 0 on  [0,  m]  and  ^ >_  0 on  [0,  it], 

but  in  this  example  it  is  not  necessary  to  do  anything  extra  to  force  this. 

Although  we  are  not  using  ordinary  algebraic  rational  functions,  we  do  have 

the  alternating  theory  in  this  situation,  and  that  is  all  that  is  required. 

To  run  this  example  we  replaced  X with  an  equally-spaced  mesh  (spacing  — — ) 

* ‘ 256 

containing  206  points.  Using  as  our  initial  reference  set  five  (roughly) 
equally  spaced  points  in  [0,  0.2"]  and  eight  (roughly)  equally  spaced  points 
in  [0.4ir,  n],  we  obtained  convergence  after  four  exchanges  and  60.0  seconds; 

||  f - r*  ||  was  1.83914  x 10  U (where  r!':  is  best).  (Note:  the  final  alternating 

set  does  have  five  points  in  [0,  0 . 2 n ] and  eight  points  in  [0.4m,  if],  but  they 
are  not  equally  spaced.)  Starting  with  eight  equally  spaced  points  in  [C,  C.2tt] 
and  five  in  [0.4tt,  tt]  , eight  exchanges  and  1 minute  11.7  seconds  were  required; 
starting  with  all  reference  points  pushed  to  the  extreme  right  of  [0.4m,  tt] 

(which  is  one  of  the  worst  possible  starting  reference  sets),  fifteen  exchanges 
and  2 minutes  0.5  seconds  were  required.  On  the  other  hand,  running  this 
problem  with  differential  correction  alone  required  5 minutes  45.5  seconds.  One 
would  expect  the  time  difference  to  increase  if  a finer  nosh  were  used. 

We  also  ran  the  Remes-Difcor  program  on  an  example  for  which  best  approx imati 
did  r.ct  exist  on  seme  reference  sets,  although  a best  approximation  did  exist 
or*  X.  Here  convergence  depended  on  the  choice  of  initial  reference  set, 
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although  we  were  able  to  obtain  convergence  ev  .n  with  a bad  initial  r* 
set  if  we  "helped  the  program  over  the  bad  spots"  by  forcing  q j>  * on  the 
reference  set;  this  (as  opposed  to  forcing  q e on  all  of  X)  did  r/.t  i *_■  qu  1 1 • 
much  additional  work. 

In  general,  the  relative  merits  of  Remes , Rtnes- dif cor , ana  difcor  for 
finite  X can  be  summarized  as  follows.  When  Remes  works,  so  does  F.o:nos-difc  or, 
and  with  comparable  speed.  Renes-difcor  will  lsually  still  work  wh  n Ren.es 
fails  due  to  problems  in  finding  a new  approximation  on  a reference  set,  and 
is  much  faster  than  difcor  if  card(X)  is  large.  Difcor  is  theoretically 
more  robust  than  Remes-difcor  since  it  does  not  require  an  alternating  tr.  :zy, 
and  || f - r^||  will  converge  to  inf  ||f  - r||  even  if  there  is  no  best 
approximation,  but  round-off  and  storage  problems  may  be  prohibitive  if 

' 

card(X)  is  too  large. 


K.. 


k initial  reference 


compute  approximation  on 
reference  sot  (Diff.  Cor. 


compute  errors,  dencniua 
maximum  error  on  entire 


We  have  t 
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Flowchart  §2.  rinding  a New  Reference  Set  (N  = n + m t 2,  NUMGR  = number  of 
grid  points) 


Let  T = tx,  ...»  tM  , t = left-most  grid  point,  t,,+1  = right-most  grid  point, 
for  i = 1,  M replace  (recursively)  t.  by  the  grid  point  between  t^_, 

and  t.  at  which  the  error  has  maximum  absolution  value  and  the  same  :'g'. 

i+l 

as  t^ , considering  only  points  which  are  not  g-poles 
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